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1  .  I  NTRODUCTION 


In  [££},  Plato  states  that  his  teacher,  Theodor  us  of  Cyrene, 
was  tht-  first  individual  to  prove  that  VS.Vfe,...  ,  Vl  7  were 
irrational  [£,  p, UoJ .  Since  then  many  scholars  have  wondered 
how  he  did  this  with  the  mathematical  tools  he  had  at  hand, 
and,  perhaps  more  puzzling,  why  he  stopped  at  the  number 

seventeen.  An  answer  which  is  almost  certainly  incorrect  but 
which  is  nonetheless  intriguing  was  given  by  J,  H.  Anderhub  in 
1918  Csee  [3,5,6,101!).  He  suggested  that  Theodor  us  might  have 
used  the  following  construction.  Form  the  triangle  with 
vertices  at  CO.Ci,  Cl, 03.  Cl, 13.  Note  that  the  hypotenuse  has 
length  Vc:.  Starting  now  at  the  point  Cl,13,  draw  a  line 

segment  of  unit  length  perpendicular  to  the  hypotenuse  and  in 
the  counterclockwise  direction;  this  gives  a  new  triangle  with 
this  line  as  one  side,  the  previous  hypotenuse  as  the  other 
side,  and  a  hypotenuse  of  length  -/S.  Continuing  in  this  way, 
a  spiral  figure  Ccailed  the  guacfratwurzelschnecAe  [103,  or 
squar e-r oot -snai 1 3  emerges  with  the  numbers  Vn  as  the  lengths 
of  the  radii  CFig.  13.  The  last  hypotenuse  drawn  before  the 
figure  overlaps  itself  has  length  "/IT,  and  this,  it  was 
suggested,  was  why  Theodorus  stopped  at  seventeen  Csee  [533. 

Although  this  is  almost  certainly  the  wrong  answer,  the 
sequence  of  points  formed  by  this  process  was  studied  on  its 

owTi  merits  by  Hlawka  CCIOD;  see  also  [£43  3.  He  proved  a 

number  cf  results  concerning  this  construction.  The  study  was 
then  pirKevd  up  by  Davis  [33,  who  noted  that  the  points  could 


Cl  .  1  J 


i  1  ^,1 

V*  *  iS.  K-  'v^t*t|.-^  >>,  <^  * 


by  the  iteratie 


2  =2  -»  »2  /  I  2  1 

r»*l  rs  r>  r> 

Cwith  •.-^  =  1 1>  and  that  they  we;  e  asymptotically  an  Archimedean 
spiral.  He  then  interpolated  an  analytic  curve  to  these 
points,  which  he  called  the  Spiral  of  Theodor us. 

Noting  that  Cl . 1 D  could  be  considered  Euler’s  method  for  the 
ordinary  differential  equation 

z  =  z^-  |z  1 


with  i!r  I  t  st.ep  size.  Davis  suggested  as  a  generalization  of 
Cl  .  1  .J  the  I  ter  a  1 1  on 


z  =  az  +  /hz:  /  I  z  I 

fh'4-i  r»  n  n 

with  o./^eC  and  z^  a  given  complex  number.  Thii  s  iteration, 
called  the  Complex  Generalized  Theodorus  Iteration  C CG7T 5 ,  is 
studied  in  [3,15,201.  Davis  later  suggested  t.he  further 
gener  al i zati on 

V  =  A«V  +  B»<V  xii''  I.  Cl  . 

n'*!  n  r<  n 

where  /  and  B  are  real  mxm  matrices.  V  is  a  cu  vt?n  nonzero 

o 

m.-vectcr,  and  ii  ■  H  is  the  Euclidean  vector  norm.  This 
iteration,  named  the  Generalize-d  Thieod.orus  Iteranor.  Py 

the  present  author,  displays  a  great  variety  of  strange 
attractors  for  appropriate  choices  of  A  and  B  CFigs.  E-5j  ,  and 
includes,  the  complex  GT3  as  a  special  case;  it  is  studied  in 
[4,15,13,191.  In  this  report  will  summarize  some  of  the 
informa  ion  contained  in  [163  and  Cl  91  C addi ng  some  material 
from  [153!1  with  respect  to  sufficient  boundedness  conditions 
and  asymptotic  estimates  for  this  iteration.  Since  w<=-  are 
inter ec  ed  in  boundedness  and  asymptotics,  in  stating  our 


results  we  will  ignore  the  case  where  V  is  such  that  V  is 

O 

the  null  vector  for  some  n  Cin  which  case  the  iteration 
stops  1>  . 


2.  Boundedness 

Let  us  try  to  find  sufficient  conditions  for  boundedness  of 
the  GTI  .  Taking  norms  on  both  sides  of  Cl  .  gives 

II V  I!  <  tl  A  li  •  II V  II  +  II B  i! 

T>*i.  r» 

Iterating  this  relation,  we  find  that 

n-  1 

11 V  II  <  li  A  i:  •  II V  II  +  II B II  •  r  II A II '  C  2 .  1  D 

n  O 

i  r  O 

and  clearly  if  llAllil  then  in  the  limit  the  norm  is  bounded  by 

IIBII/Cl  -  IIAIO 

Hence  iAII<l  is  sufficient  for  boundedness  of  the  orbit  C  f  or 
any  •  From  this  and  the  similarity  of  Cl  .  £7^  to  the  linear 
syst-em 

V  =  AiiV  C~  2? 

n-^l  r> 

which  fas  a  bounded  solution  for  every  when  the  spectral 

radius  of  A  is  less  than  unity  Cand  at  least  one  unbounded 
solulicr.  when  the  spectral  radius  of  A  is  greater  than  uni  tyT;  , 
it  seems  reasonable  to  conjecture  that  pCATXl  suffices  for  the 
boundedness  of  Cl.  27),  where  pC  AT)  is  the  spectral  radius  of  A. 
That  tf.is  is  in  fact  the  case  follows  from  two  lemmas.  We 

omit  th >  rather  technical  proof  of  Lemma  1.  It  is  not  hard  to 
establi'.h  this  result  by  using  an  approach  similar  to  that 
ih  live  SkPcI  six  but.  it-  is  rather  tedious  to 

carry  c  it  the  F>roof  in  detail  Csee  [15]  7).  Lemma  2  is  a 
discrete-  version  of  an  ordinary  differential  equations  result 


4 


^  .w.,  ^  r 

Lemma  1  •  If  li  V  li  i  5?  unbounded  -then  II V  1*' — >o:>. 

n  r), 

Lemma  2  Let  A  b»  a  matrix  with  pC AD<1  and  suppose  that  BC n j 
is  a  sequence  of  matrices  with  HBCnjlf — »0  as  n — *oo.  Than  ali 
solutioriS  of 

z  =  C  A  BCnJJ  -z 

n*l  n 

tend  to  zero  as  n — »<xi. 

Proof:  Since  pCA!)<l  there  exists  a  matrix  norm  li  •  II  such  that 

A 

HAH  <1.  We  will  choose  this  to  be  a  natural  matrix  norm  and 

A 

Will  a]  SO  denote  by  I;  •  li  the  vector  norm  which  induces  it. 

A 

Lee  e  equal  Cl-liAli  J .  Now  since  all  matrix  norms  are 

A 

equivalent  in  the  sense  that  if  liM-.lnl)ll  — >0  then  liM'ntN  — >0,  we 

o  ft 

have  from  our  hypothesis  that  liBCnDH^ — *0.  Hence  we  can  choose 

an  N  such  thac  i<  B<.  i  j  li  %  e'*‘c:  i  or  al  1  i  cCN. 

A 

Now  consider  iterating 


z  =  CA  +  BCn'i'iT 

r> 

C  f  or  sore  z  li  N  times  to  produce  the  vector  z  .  Set.  o'  =z  anr 

O  NON 

f''  —  C  A  +  N  +  nii  '■>  • 

r>*l 

Vnl'O.  Clearly  /.*’  =z  for  all  Now 

I  N*\ 


o>  -  C  A.  BC  N  +  nJj*  ...  *■;  a  +  BCNll)*"' 

O 


and  so 


li  w  II  <  II A  +  EC  N  +  rft  li  « 

n+J  A  A 

rs 


*iiA  +  BCN5II  *11  u>  li 

A  O  A 


=  li  W 


1  II  n  H  A  +  BC  N  +  i  j  li 

O  A  '  '  i 


men 


tl  CO 


A 


<  llw  !l  n  ^  *IAil 

O  A  *  *  A 

V  =0 

r> 

<  ii  lo  11  n  ^  1 

O  A  *  * 


I  =o 


lIBCN  i:>(I  3 

A 


:=  11(0  II 

O  A 


n 


wHi  rh  -leariy  tends  to  zero  as  n — »flD.  Hence  llo  il  — »0  as  n — »a' 

n  A 

so  the'.  l!z  ii  — >0  as  n — kxi,  and  so  by  the  equivalence  of  vector 

»h  A 

norms  on  [R"’  we  have  that  Hz  il — ►©  as  n — koo.  This  completes  the 

Th 

pr  oof .  D 

We  may  now  show  the  desired  result. 

Theorem  1:  If  pCAri<l  then  the  solution  of  C 1 . 23  is  bounded. 

Proof:  Suppose  the  contrary  and  choose  a  V  such  that  II V  H  is 

O  n 

unbounded.  Py  Lemma  1,  11 V  II — *et>.  Define  BCn3=B/IIV  ii .  By  the 

Th  n 

hypoth^s s  ,  IIFCn3ll — >0.  Consider  the  iteration 

z  =  CA  m  BCn3  3*z 

rv*1  n 

With  z  =V  .  Clearly  z  =V  for  all  n.  The  conditions  of  Lemma 
<>  o  ^  r> 

2  are  met  and  so  llV  II — >0  as  n — »oo.  But  this  contradicts  IIV  Ii 

r>  n 

unbounded,  which  establishes  the  result,  o 

It  also  possible  to  show  that  if  pCA*>>l  then  there  exists 
a  V  s'.jch  tl'.at  nV  II  is  unbounded  Csee  section  53.  In  ohis 

O  n 

sense  'he  GTT  is  sur  pm  i  si  ngi  y  similar  to  the  linear  system. 
C  2 .  23  . 


In  tJ  .= 

part.icular  case  A=v3  a 

del.ai  1  ed 

anal ysi s  is 

C 16. 1 7]  , 

Th):  s  case  is  relat.ed 

to  the 

power  method 

for  t-he 

numerical  solution  of  eigenvalue  problems. 

3.  Asymptotics 

This  leaves  the  question  of  boundedness  for  the  case  pCA3=l. 
It  IS  possible  to  construct  examples  showing  that  when  A  is 
uni  spec '.ral  the  orbits  of  the  GTI  may  be  always  bounded  for 
one  chc.,  ce  of  b  and  always  unbounded  for  another  choice  of  P. 
However,  a  general  statement  about  the  asymptotics  c^f  the 


(=. 


4  ♦  •♦•lii. 

j  I-'*-  7 


*  i  t:**'  ***'s**  y'r:^  U'Siii*-'  A  i  »r  V,/-’,»  4  ♦•♦•lii.  ^  ^  r,  *  i  »r 

^  S.-  V»  i  CA  V‘  >_<’  J  t  J"  •w  i4l<^  S.AV.-  *■  A  i  ^  V-  ^  7  ^«'  i  *  CA  S.  aLw^7 

pCAD<i  and  any  eigenvalues  of  A  that  lie  on  the  unit  circle 
are  simple  Ci.e,  have  as  many  linearly  independent 
eigenvectors  as  their  multiplicity;).  In  this  case  the  linear 
system  C2.2D  has  the  property  that  all  solutions  are  bounded 
Cal  though  in  general  the  bound  will  depend  on  ^  ■  Before 
statin-c  the  corresponding  result  f c-.>-  the  GTI  ,  .'fe  state: 
Hukuwara's  Theorem:  Let  be  any  solution  of 


:  =  CA  •+  Brt5:n«z 

t*i  t 


C3.  1  D 


for  t=C',l.c;..  .  .  where  A  js  a  constant  matrix  and 

CT' 

J]  liPCtJll  <  00 
t  =  1 

Suppose  that.  A  is  of  bounded  type,  so  that  all  solutions  of 
y  ~A*y  are  bounded  for  t>0.  Then  every  solution  of  C3.i;) 

t  '  i-i 

IS  hounded. 


Thii  s  result,  an  analogue  of  a  theorem  f  romi  differential 
equations:,  may  be  found  in  [211.  The  proof  there  uses  the 
v«=-ctor  norm,  but  the  argument  is  valid  for  any  vector  nor  rri  and 
compatible  matrix  norm.  We  will  now  use  it  to  prove  an 
asymptr  1. 1  r:  result  for  the  GTI, 

Theorem  2;  Let  A.  be  of  bounded  type  m  C 1 , 27^ .  Thien 

OD 

H  IXliV  II  =  CD 

r> 

7  =0 


Proof:  Suppose  the  cont’^ary,  i  .  e.  that  for  some  V  we  have 

00 

J’,  1  /II V  II  <  CD  C3.  21) 

r> 

n  =  0 

and  construct  such  a  sequence  <V  >  by  using  this  V  and  Cl  .  27' . 

T>  O 


Define  3Cn7=B/llV  11  and  consider  the  iteration 

n 

Y  =  CA  +  BCnjDRY 
n-t-l  n 

where  ve  set  Y  =V  .  Clearly  1'  =V  for  all  nLO;  the  two 

00  r>  rt 


Now 


sequences  are  identical. 

00 

Y.  liBCn:)  II 

r>  =  O 


00 

Y  liBII/ilV  II 

r\ 

=  O 

CO 

tiEi;  Y 

r> 

n  =  O 


which  by  supposition  is 
Theorer,,  HY  li  is  bounded 

r^ 

such  that  IIV  IKM  for  all 


bounded.  Therefore  by  Hukuwara’s 
Then  IIV  11  is  bounded,  that  is,  3M>0 

r> 

n.  But  then  111/^/  ll>l/M  and  since 


n 

l/llV  II— x-»0  we  cannot  have  C3.C.3. 


n 

Thiis  contradiction 


n 

establishes  the  result,  a 

ill  fa.ct  a  slightly  more  general  result,  holds.  If  It  •  li  is  any 
vector  norm  and  5C  •  iJ  is  any  function  that  is  continuous  on 
rO.arO  and  positive  on  CO,od:;,  and  if  A  is  of  bounded  type,  then 
any  sol  ut. ion  o: 


V  =  A*v  +  E*V  /  JCliV  ir^ 

n-^l  n  n  n 

must  be  such  t-hat 

cr. 

Y  1  II V  IID  =00 

r\ 

r>  =  O 

Thiis  mta  •/  be  shown  in  a  completely  analogous  ma.nner  C  see  Cl  53 


f  or  det.^ils  and  appl  i  ca ti onsl  . 

TbiC-orem,  £  establishes  that,  when  pC  A.j  =1  but  A  is  of  bounded 
type,  if  the  orbits  diverge  then  they  do  so  in  a  way  that  is 
not.  much  worse  than  linear  C  si  nee  the  sum  of  the  reciprocals 
of  the  normc  diverges!).  In  fact  if  ilAll=l  in  the  spectral  norm, 
then  it.,  is  a  fact  that  the  divergence  is  at  most  linear,  as 
can  be  seen  by  setting  II A  11=1  in  C2.  1>  to  get 

ilV  It  <  IIV  It  +  nllBlI 

r»  O 

as  the  asymptotic  estimate.  We  will  later  consider  this  point 
in  grea  er  detail  when  we  show  in  secti  on  5  that  if  A  i  s  of 
bounded  type  then  the  divergence  is  always,  at  most  linear. 


>_•.*  tTvtl  ^  i  IN'-  Cr_^Oj*N.-civ  4.  •-•»i 


i  n  or  i-i  ’—r 

refine  our 

as  ym.pt ot  i  c 

estimates  of  the 

previ  ous 

sect  1  or , 

we 

wi  i  1  need 

a  for  mul a 

for  the  solution 

of  the 

1 terat j  on 

C  1 

.  £3  .  Gi  von 

an  1  tor at ion 

1  of  fho  f  or  m 

2  =  ACtJ^^Z  +  Cxi 

t  t-l  1-1 

With  c  a:  d  2  Qi  v<^r*  vectors  and  AC  1 1  a  time-varving  matrix, 
V  o  ~ 


ttic”  5‘-_'iWtic*r,  i._ari  a^-wa  3'*^  htr  wr  itt-en  in  thf  2  orm 

t 

z  =  r  Y  -  Y'^  o.'  +  Y  •y"‘-c 

t-**!  IS  S  O 

s  =  C» 

wTiere  c  ^  t.he  initial  condition  2  ,  and  -CY  >  is  a  f  undamentai 

o  t 


m=trrax  set  of  solutions  for 


z  ~  As.  t  j 

i  t-l 


Qh.i") 


Cthat  IS,  Y  =I  and  Y  =ACnl>*Y  D  [  £i  3  .  For  the  GT7  ,  the 

O  r»*l  n 

matrix  A  is  const  snl  ,  a  no  z  =V  is  given.  Furthermore,  Y  =a'" 

O  O  r, 

for  the  linear  system  C£.c:t  corresponding  to  C4,1Y,  and 

f.>  =  FxV  X  I.  V  ii 

t,  t  t 

SO  that 


a’^**V  +  T  At  *«BxV  /  |iV  |i 

o  ^  t  t 


C  4  .  aY 


is  the  solution  of  Cl. 13.  Tni  s  is  called  the  Green’s  Function 
Representation  of  the  solution.  We  will  use  this  form  of  the 
GTl  to  address  the  questions  of  boundedness  and  asymptotics. 


5.  A  Direct  Approach  to  Boundedness 


Thie  apf  'oach  used  in  the  se  :ond  section  to  show  that  p<:h'><l 
suf  facet  for  the  boundedness  of  C 1  .  £3  is  trie  one  coiginaily 


ompl  oywd  by  author  in  f  j  5j  and  also  appears  in  [lb],  A 
more  d-.rect  method,  however,  was  used  by  the  author  in  flQ], 
and  we  now  give  this  method 

Suppose  that  pCA3<l.  Then  there  exists  some  natural  matrix 


norm  1'  • 

H  such  th 

a 

at  II  All 

a 

<  1  . 

Us  i  ng 

II  ■  il 

a 

to 

r  f=9pr  esf^nl 

the 

vector 

norm  which 

i nduces 

thi  ^ 

ma  t  r  i  X 

nor  m 

as 

w©l 1  ,  we 

have 

from  C4 

.  Si  that 

IIV  li 

n+l  O 

<  H  a’"**V 

o 

II  + 
a 

1  E  A-' 

*«BxV  X 

i 

i:V  |! 

1 

i 

i 

id 

t  =o 


VndO.  Thus 

rv 

iiV  li  <  II  Ai: •  liv  I!  +  r  llAl^'-llEli  -i. 
n+j  C*  0(  0  0  0(01 

t  =o 

where  t  is  a  positive  constant  such  that 

II Vi;  xliVli  <  k 
o 

for  all  nonzero  vectors  V.  Such  a  k  exists  since  all  vector 
norms  are  equivalent  on  CR”’  CISj.  Summing  the  geometric  series 
and  col  ler  ting  like  terms  gives 


iiV  |i  <  <eC  ni  + 

k  •  ilBi 

XCI  -  llAli  1 

>-i*l  o 

Cx  0« 

wner  t  c-. 

.i-.  ts.!  as  n — Hence 

HV  II 

IS  asymptotical] 

y  bounced , 

ry  O 

and  by 

the  equivalence  of 

vector 

norms  so  is 

li  V  li .  An 

r* 

asympi  o- 

ic  bound  for  H V  II  is 

r*  O 

g  i  ven 

by 

k-llBl!  XCl  -  HAH  j 
c* 

in  t.errrs  of  the  o-norm.  When  the  spectral  norm  of  A  is  less 
than  unity  the  corresponding  bound  using  that  norm, 

ilBlIxCl  -  II  AID, 
is  gene -ally  fairly  tight. 

We  have  shown  that  pC  Al)  <  1  suffices  for  boundedness.  In  an 
anal  ago  IS  way  we  can  show  that  if  p(TAr>>l  then  there  must  exist 
some  V  such  that  |iV  il  is  unbounded. 

1  0 


C  A  s 1  mi 1 ar  appr  oa c  h  has 


properties  similar  to  those  of  the  linear  system  C2.2I>.  In 
the  next  section  we  will  derive  improved  asymptotics  for  the 
i  ndeter  ;ai  nate  case  where  A  has  unit  spectral  radius.  Before 
doing  Sf'  lei  us  collect  the  information  above  as  a  theorem. 
Theorem  3:  Consider  the  iteration  C 1 . 25 .  If  pCATii'l  then  the 

iteraticn  is  bounded  for  any  V  ,  and  if  pCA5>l  then  there 


11 


®xa,  stfc 


chat  the  iteration  as  unbounded. 


5.  Improved  Asymptotic  Estimates 


From  C5.  it  is  immediately  ap<parent  that  when  the 

divergence  will,  in  general,  be  exponential.  This  again 
leaves  the  case  pCAj=l.  As  before,  let  us  first  look  at  the 
case  where  A  i  s  of  bounded  type.  Using  the  spectral  matrix 
nor  m ,  C  4 .  5  j  gi  ves 

i;  _ _  i; 

C6.  1  ^ 


II V  u  <  |iA""‘v  li  i  r  ‘xBxv'  /i;  V  li  " 
r\*i  O  II  ^  t  t 

1=0 


<  K'  -t-  E  IIA 
t  =o 


IBi. 


where  K  2:0  is  a  bound  on  HA'  *V  II  C  which  in  general  will 

i  O' 

depend  on  .  Now  since  A  is  of  bounded  type,  liA^II  is  also 
boundec .  that  is,  IIa'^II<K  for  some  K>0.  Thus 


II  <  K  +  IIBII  r  K 

►i  1  ^ 


t  =  o 
<  K  +  K  n 

i  z 

showinc  that  the  divergence  is  no  wc-rse  than  linear. 

If  pC  A'J  =1  but  A  is  not  of  bounded  type  then  it  is  possible 
to  show  that  li A*^ ii  =€>C n"'D  Cin  fact  IIA'^ll=©Cn^  where  p<m  is  the 
maximal  degree  of  all  nonlinear  divisors  of  A  associated  with 
unimodular  eigenvalues!),  where  A  is  mxm.  Hence  from  C6.17i  we 
have 

IBIi 


liv  II  <  II  a'"**  II  •  II V  II  +  r  iia’i 

n+1  O  ^ 

i  -O 


and  clearly  II V  li  diverges  at  most  as  CJn  fact,  it 

diverge^:  as  OC  ,  and  hence  in  the  worst  case  as  OCn^^^. 
Hence,  when  pC  A”^  =1  ,  if  the  iteration  diverges  then  i  i-  does  so 


1  r- 

X 


',■«-•  li. '-•<-;  j  I  f.'ti'*  A  j  .  »»‘.  •*  »  A  •.- -s.-  K'»t^  w  o.  _'  o.  >•-  i.i 


Theorem  4; 

Uonsi der  C  1  . 

03  and  suppose  that 

pc  A3  =1  and 

t  he 

i ter  at 1  on  i s 

unbounded . 

If  all  eigenvalues 

of  A  wit  h 

uni  1 

modulus  are 

Simple  then 

the  divergence  is 

no  worse 

t  han 

linear.  C>t.herwise  the  divergence  is  no  worse  than  polynomial 
of  order  m,  where  mis  the  order  of  A. 

7.  Li  mi  t  Sef  s 

Let.  !.is  ri^f.ine  TCA.B.V  to  be  the  limit-  set  .'.or  c--i  i  m;  t  set 

o 

[15]3>  of  the  iteration;  that  is,  YerCA.E.V  1>  if  there  exists 

o 

an  1  ncs (=•=.<:  1  n.-!  subseauence  nC  1 1;  such  that  V  — »Y  as  i — ►a'  when 

'  r»(  t  > 

the  i  n.t  t  1  a  j  condi  t.i  on  is  V  ,  Two  t  heorems  f  oi  ]  ow  eas  .  1  y  f  r  cm 

(I 

t.  he  ma  t  =  r  i  a  .1  in  [133. 

Theorem  Z\  St  ipfpos^  vrjds*,.  1  or  sojTic"  \  1  sscj^terriCt:  "v.  v  ^  is 

O  rt 

bounded  Thien  FCA.B.V  D  is  a  rion-eiriDt  v  comoact  set  and 

o  ... 

V  — »rCA,B,V  r>  as  n — »oo. 

Ti  O 

Broof;  til  ear  jy  the  closure  of  P  is  simply  P  and  so  P  is 
closed.  Tne  bounded.ness  of  fV  >  implies  that  P  is  also 

r> 

bounded  Kence  P  is  compact  C  by  def  i  ni  ti  ori[<  .  By  the 

Boi zanc -Wei erstr ass  Theorem  [233,  P  is  not  empty. 

Now  V  — »P  means  that  infCliV  -Y II :  YeP  j — >0  as  n — mX’  [133. 

rs  rv 

Suppose  the  contrary.  Then  we  can  find  a  subsequence  of  <V^> , 
indexed  by,  say,  nCi3,  which  remains  a  finite  distance  from  P. 
Now  V  is  a  bounded  infinite  sequence  and  so  by  the 

r- 

Bol  zano-Wei  erst  r  ass  Theorem  it  contains  a  limit  point, 
that  pc  nt  must  be  in  P,  which  is  a  contradiction,  n 

1  5 


But 


Theorem  C.,-  Suppose  that  for  some  V  the  sequence  iV  >  js 

O  r. 

houfide  o ;  ,c  0<?'~C  A  ,  F- .  \  'j  .  Then  PC  A.  F' ,  V  'i  is  po<--  i  l  i  vel  v 

o  o' 

i  nvar i ant . 

Proof;  Def i ne 

TCVO  =  A«V  B«V/(iVli  C7.1D 


so  that 


V  =  TCV  D 

n-*-!  n 

is  the  iteration  under  consi derati on .  Since  OeT.  the 

t r ansf c r ma ti on  T  is  continuous  in  some  neighborhood  of  any 
point  Ver.  Fix  Y  and  choose  a  subsequence  nC  i  ■>  such  that 
V  — »Y .  By  the  continuity  of  T,  we  have  that. 

TCV  :)—^TCyj 

tM\t 

as  i — »iT  But 

TCV  j  =  V 

rd'  r<d>*l 

SO  that 

V  — ♦TCYD 

wher  e  mC  iJ=nCi!)-*-l  .  Hence  TC  YD  eF ,  a  riW  r  i  3  pos  i  1 1  vel  v 

invariant,  d 

In  fa  1..  It  IS  clear  frnrr.  t-hie  above  that  TtyteF  whenever  YeF 
and  TCY!)  is  defined.  In  other  words.  rxco>  is  positively 
invariant.  Further  topological  results  on  discrete  dynamical 
systems  that,  can  be  extended  to  this  ca'^e  Cwith  some  care 
taken  near  the  origin!)  may  be  found  in  [141. 

Suppose  that  YeFCA,B,V'  j,  so  that  V  — ►Y  for  some 

O  rid) 

increasing  subsequence  nCiD.  Since  the  transformation 

T;  [F  — >[F  of  C7.ll)  has  the  property  that 

TC-VI)  =  -TCVj 
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t  ha 


o 


ruv.  > 


it  fellows  that  the  sequence  qenerated  by  -V  is  -V 

tends  ■  -n  -Y  ■  that  is.  -Y<=rrA.B.-V  '> .  In  a  similar  w-jin 

o 


v-li.  c!. 


easv^  t  c.  sec  that  i  ei.  A*  £:>,  y  anq  ^  tn.  i  A.  iii.  V  ^  . 

o  o 

Xr.  tr.e  case  ol  a  qlobal  attract  or  •  F'cA.&.V  j  is  xr.dleoer 

'  o 

ol  *  1  or  the  oar  oi  cui  ar  A  anq  q  unoer  consi  der  aci  or..  ^  r. 

o 


C  , 

i  t 

foil ows 

f  r  om 

the  above  that  the 

ClC^t,rOP 

central 

i  y 

symrr.etr  i  c 

.  1  .  e  . 

YeF  implies  -YeF.  and 

that  the 

attract 

OP 

is  f  ound 

when 

A  and  E  are  replaced 

by  -A  and 

respec i  vel  y .  Some  of  the  attractors  displayed  bv 

iteration  f  r-r  various  choices  of  A  and  E  are  shown  in 
iiqurei  i  see  aiSO  i^j.  iser^ies  has  found  ar.  anaxVt 

expression  for  the  attractors  for  a  certain  set  of  case 
the  lorn.  A=or.. ,  o'lC  Csee  C411/.  For  the  case  A  =  0  tr.e  for 
the  attractor  is  also  known  C15.15,1F‘. 


1  s 


same 


-B. 
tht  s 
the 


Ol 


m.  of 


on 


The  Generalised  Theodores  Iteration  is  in  some  wavi  very 
siiTiJxa.-  t. t-i"ie  sijT.q.ie  iinecP  system,  i  xt"iec*rem'  ar.o  ^et  it 

contains  some  very  complicated  structures  CFigs.  £-51o  It  is 
likely  that  the  attractors  are,  in  general,  strange  but  net 
always  chaotic  Csee  [9]3>;  some  computational  evidence  to  this 
effect  has  been  provided  by  James  Heyman  of  the  Naval 
Postgraduate  School  in  his  study  of  the  use  of  chaotic 
discrete  dynam.ical  systems,  particularly  the  Henon  attractor, 
as  pse  idor ar.dom  number  GS'r’ifcr'r  o  X*  or  s  C  soos.  1  c  w^rnTr>^^r*ii  c  ^  i  c-r"*  v*  . 


i  W  -  w  ^  1  w  . 

I  v..*iC  ^  «  A  Vw  * 


oi  tr  X  Tk  or*  Odr  • 


disf  i  ni  1. 1  ons  of  the  terms  "chaos"  and  "fractal"  Csee  [  1  1  it 

IS  difficult  to  prove  anything  substantati ve  about  the  chaotic 
r  i  a  t  u  r  a  c^l  t  ti  a  a  1 1  r  ol  o  t  *  r  s  b  u  <*-  >  d  s  n  tw' t  in  rna  n  ^ 

w’el  1  K'novn  strange  attractors  ara  still  awaiting  artal  i  c  a  ^ 
verification  of  the  observed  qualitative  features  of  the 
underlying  iterations  and  differential  systems.  Orie  aspect  of 
the  GTT  which  wo  feel  is  particularly  relevant  to  the  field  is 
the  v'a  s  t  a  r  r  a  y  ot  geom.e  tricai..ty  distin^-i  attract.  *-’rs  i-o  be 
founo  i.n  the  lor  dj.>ierent.  chc>ices  ol  the-  miatrices  A  and 

E;  even  thougr.  these  two  matrices  essentially  represent  eight 
parameters.  the  number  of  distinct  attractors  is  still 
s  ur  pr  i  1  r*o  ♦  c*.  nc*  o-r  i?  nc*  w  &  w^r  o*  c».r4C>l*»*rr  ^  P-,«r*or 


m.ap  With  this  property. 


Further  i  nvesi.i  pat  i  on  of  this 


p's*  v..  '  j  k  J  O 


seems  warranted. 


Hopeful  1 y  i he  f i cur  es  c i ted  her  e 


and  in  ilEl  will  miotivate  research  on  tl'.e  question  of  the 
number  of  distinct  attractors  which  can  belong  to  a  single 
j  te  r  at  1  rn.  as  t  he  'worl  s  bv  Gieicl  C  ^  ^  .  Eevane^^  i.  s-..-  .  and 


Iftigi^ins  [ciOfCiUj  have  Ser  ve..i  tv.*  hei  p*  m*..>  v  ^  va  t  e 


author . 
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